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ABSTRACT 

We examine the linear density contrast at collapse time, <5 C , for large-scale structure 
in dynamical dark energy cosmologies, including models with early dark energy. Con- 
trary to previous results, we find that as long as dark energy is homogeneous on small 
scales, § c is insensitive to dark energy properties for parameter values fitting current 
data, including the case of early dark energy. This is significant since using the correct 
5 C is crucial for accurate Press-Schechter prediction of the halo mass function. Previ- 
ous results have found an apparent failing of the extended Press-Schechter approach 
(Sheth-Tormen) for early dark energy. Our calculations demonstrate that with the 
correct 5 C the accuracy of this approach is restored. We discuss the significance of this 
result for the halo mass function and examine what dark energy physics would be 
needed to cause significant change in S Cl and the observational signatures this would 
leave. 
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1 INTRODUCTION 

Observations of Supernovac 
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la llRiess et alj 1 19981 ; 
Kowalski et all l2008l).~ the Cos- 



iPerlmutter et all [l99!_. 

mic Microwave Backgroun d (CMB)dKomatsu et al. 
and large scale structure (|Cole et aL 20051 ) demonstrate 
that the expansion of the Universe is accelerating. Dis- 
covering the physics of the dark energy thought to be 
driving this phenomenon is a key goal of modern cosmology. 
Observations that probe the non-linear growth of structure 
are sensitive to the entire history of the Universe and are 
a crucial element in attempts to measure the evolution of 
dark energy properties with time. 

The abundance of collapsed structures as a function 
of mass, the Halo Mass Function (HMF), is an impor- 
tant stati stic that is measurable through strong lensing 
statistics ( Bartclma nn et al.l 1 19981), galaxy redshift su rveys 
|Evrard et all 120021 ) and X-rav l|Borgani et all 1200 ll ) de- 
tection of clusters and future cluster surveys utili sing the 
Suny aev-Zel'dovich effect signature in the CMB (iTauberl 
2005). Accurate estimation of this statistic as a function 
of cosmology is therefore required in order to extract the 
maximum amount of information from observations. 

Determining the HMF of cosmological models can be 
a computationally expensive task requiring many large N- 
body simulations, since non-linear gravitational structure 
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growth cannot be calculated analytically. Therefore, sim- 
ulation calibrated tools for rapidly and accurately generat- 
ing the observational signatures of a wide variety of models 
are essential. However, given the theoretical uncertainty sur- 
rounding the physics of dark energy, tools for predicting the 
HMF must be valid for dark energy models models as gen- 
erally as possible, to avoid detailed computation for every 
one of the plethora of possibilities. 

Current method s for estimating the HM F fall into two 
main categories fsee ICoorav fc Shcth (2002 ) for a r e view) . 
The first are methods based on Press fc Schechterl (Il974l ) 
theory that relate the density of collapsed objects of a given 
mass to variance of the density field on scales enclosing that 
mass in the mean, a 2 (M), and a threshold parameter, S c , 
determining the linear overdensity required for collapse by a 
given redshift. The lead i ng ap proach based on these ideas 
is the ISheth fc Tormenl l| 19991 ) (hereafter ST) mass func- 
tion, which incorporates ellipsoidal collapse (rather than the 
purely spherical collapse of Press-Schechter theory) and has 
free parameters that are simulation calibrated. 

The second type of mass function fitting approach is to 
fit directly the multiplicity function 
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in which P(k) is the power spectrum of density fluctuations, 
p is the mean matter density, and W(k, M) is the Fourier 
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transform of a spherical top hat function with a radius that 
encloses the mass M at the mean density of the universe. 
Ijenkins et al.l (|200lh (hereafter J01) found a universal func- 
tion of a from simulations of ACDM, matter dominated and 
open universes that fits the multi plicity functi o n for all of 
these cosmologies. More recently IWarren et al.l l|2006l ) have 
found a similar formula that fitted their simulation results 
slightly better. 

Since there are many different dark energy models 
currently proposed, it has not proved practical to include 
dark energy variation into HMF fitting work, and cur- 
rently available formulas are calibrated only to ACDM and 
matter only cosmologies. However, the J01 formula has 
been demonstrated via N-body simulations to be valid (at 
least at the 20% level) for an evolving dark energy model 
(ILinder fc Jenkins|[200a ) and the ST formula was also found 
to agr ee with evolving dark energy models in iKlypin et alj 
(2003), indicating that the universality of gravitational col- 
lapse appears to extend to at least some dark energy models. 

One particularly interesting class of evolving dark en- 
ergy models is the conc ept of early dark energy (hereafter 
EDE ). In these models (|Wetterichll2004l ; iDoran fc Robbers! 
2006), dark energy has a non-negligible energy fraction 
through the entire course of cosmic history, rather than 
being important only at late times as in the ACDM 
model. In an investigation of the HMF in EDE cosmology, 
iBartelmann. Doran fc Wetterichl (|2006l ) (hereafter BDW) 
calculated, using the ST approach, that the presence of EDE 
leads to a significant enhancement of the abundance of col- 
lapsed objects relative to ACDM, particularly at high red- 
shift. However, in a recent study l|Francis et al.ll2008l ). we 
found that the J01 and Warren et al. mass functions, in 
contrast, predicted much less difference between EDE and 
ACDM and N-body simulations agreed with this result. The 
ST results using the BDW <5 C were clearly ruled out by sim- 
ulation data. 

Why did the Press- Schecter (Sheth-Tormen) approach 
appear to fail? While the multiplicity function is expressed 
directly in terms of a, the Press-Schechter approach involves 
5 c /a, suggesting that 5 C could be the source of the discrep- 
ancy. BDW calculated that S c is significantly altered in EDE 
cosmology compared to ACDM and this difference in S c re- 
sulted in the difference in mass functions. iGrossi fc Springel 
(2008) noted that if instead of using spherical collapse ar- 
guments, a constant value of S c = 1.689 is assumed for all 
cosmologies at all redshifts the basic agreement between the 
ST and J01 mass functions is restored. In this work we inves- 
tigate this issue further, re-examining the root calculation 
of 5 C . 

In Section [2] we re-examine methods for calculating S c 
and advocate a more accurate technique for dark energy 
models, especially for early dark energy. Given our calcula- 
tions of 8c, in Section [3] we determine halo mass functions 
from the ST approach and compare to the J01 approach. 
We also discuss the implications for mass functions in dark 
energy cosmologies more generally. Finally in Section [4] we 
discuss the dependence of 5 C on dark energy and suggest 
under what conditions we might expect a significantly dif- 
ferent value from that of ACDM and the observational con- 
sequences. 



2 COMPUTING THE LINEAR DENSITY 
CONTRAST 

The picture of non-linear growth of structure involves den- 
sity perturbations growing in amplitude, initially by a linear 
growth factor, then achieving sufficient density contrast to 
separate from the Hubble expansion and to collapse, increas- 
ing its density in a non-linear manner. One can calculate 
the level to which the density contrast would have grown 
in linear theory by collapse time, as a convenient parameter 
(and an essential ingredient in Press-Schechter formalism), 
though the true, non-linear density contrast is much larger. 
The linear density contrast at collapse, S c , is defined by 



« e =]im[A(a)-l]£^ 



(3) 



where A is the overdensity, p/p, of some spherical region 
of the universe that collapses at scale factor a c , D+(a) is 
the linear growth factor. This parameter, therefore, quan- 
tifies the linear growth from the early universe until a c of 
an overdensity known to collapse under non-linear growth 
by that time. For matter dominated cosmologies, S c = 1.686 
and for ACDM it becomes a weak function of cosmology, 
retaining the matter dominated value at high redshift and 
dipping only slightly by redshift zero. 

In order to calculate <5 C , we need to solve the linear 
growth equations to obtain D+(a c )/D+(a). We also need 
to find the overdensity in the early universe that leads to 
collapse of the perturbation at the e xact desi r ed sc ale fac- 
tor a c . The usual approach (see e.g. IPeeblesI fl980)) is to 
solve simultaneously the Friedmann equations for the back- 
ground universe and perturbation, treating the perturbation 
as a closed universe. It is also traditional to normalise these 
equations to turnaround, when the time derivative of the 
perturbation radius is zero. Following the notation of BDW, 
we therefore need to solve the following equations (for sim- 
plicity we take a flat background universe) 



x — \j h \x 2 g(x) 

x 



y = 



ujC, 1 + 3w(x) 
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where x = a/a ta and y = R/Rt a , and R is the radius of 
the perturbation (which is assumed to be spherical). The 
dimensionless density parameters of matter and dark energy 
at turnaround are ui and A respectively and ( quantifies the 
overdensity at turnaround; w(x) is the dark energy equation 
of state and g(x) is the dark energy density normalised to 
the turnaround value. The dots indicate derivatives with 
respect to the time parameter r = Ht a t where H is the 
Hubble parameter and t is cosmic time. Finding the value 
of C, that ensures collapse of the perturbation at the required 
scale factor a c requires a numerical search. 
The overdensity, A (a;), can be defined via 



c* 3 



(6) 



and therefore, once £ is found, the behaviour of Eqs. @ and 
(J5| determines the size of the overdensity at early times. 
In BDW, an approximate solution for the overdensity 
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at early times was derived of the form 

A(x) = 1 + -F( 1/3 x 
5 



where there are two solutions for F, one for the case when 
the dark energy equation of state at early times, Wi = 
liniz^o w(x) < —1/3, and one for Wi > —1/3. ACDM is 
an example of the first case, and the solution results in 



Facdm = 1 + 



(8) 



The second case includes the IWetterichl (|2004h EDE model 
examined in BDW, with 



Fede = 1 — 



i-n e 



(9) 



where D. e is a parameter of the lWetterichl <|2004h EDE model 
quantifying the fractional dark energy density at early times. 
In this model, as Q e — > 0, the cosmology converges to 
ACDM, and hence we should expect that these two solu- 
tions also converge. However this does not occur, in this 
limit Eq. © converges to unity and not to Eq. @. 

Moreover, the non-negligible presence of dark energy 
in the early universe should, relative to ACDM, slow the 
collapse of the perturbation due to the lower matter clus- 
tering source term in the early universe and later the higher 
expansion rate. We should expect then, that for the same 
collapse scale factor a c , the addition of EDE should require 
increasing the overdensity in the early universe in order to 
compensate for the slower non-linear growth. However, the 
solution for A(x) in the early universe from BDW predicts 
a decrease in the overdensity for EDE relative to ACDM. 
These problems raise doubts about the accuracy of this so- 
lution. 

We note that a precise determination of the initial A(aj) 
is needed in order to accurately determine S c . We find typical 
values of A(aj) to be A(ai) ~ 1 + 3(ai/o c ) where, to ensure 
numerical convergence (results independent of the chosen 
aj) we take < a c x 10 -4 . This means that even a small 
error, say of order 10 -3 in A(ai), is greatly amplified when 
going to A(dj) — 1, needed to calculate S c in Eq. Q. In 
this example the error induced in S c would then be of order 
unity. 

We employ a different, straightforward approach to 
computing S c , purely numerically. We still solve Eqs. ((4]) 
and ©, however rather than searching for the overdensity 
at turnaround numerically then scaling this back to early 
times approximately, we perform the numerical search start- 
ing the integration at early times, and hence search directly 
for the overdensity A(di) that causes collapse of the pertur- 
bation (y — > 0) by the desired collapse scale factor a c . Our 
approach has several advantages: 

• In the previous approach, a numerical search must be 
made to find the overdensity at turnaround before this can 
be scaled back to the early universe using an approximate 
solution. Performing this search at the early time instead 
avoids the need for the approximate scaling back without 
adding to computation time. 

• Finding the overdensity at turnaround implicitly as- 
sumes that the rise and fall times of the radius of the per- 
turbation are equal. This is not true in general dark energy 
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Wetterich (2004) EDE Model numerical solution 

Wetterich (2004) EDE Model BDW solution 

Doran & Robbers (2006) EDE Model numerical solution 



Figure 1. Linear density contrast at collapse time, 5 C; for ACDM 
and EDE cosmologies. The ACDM model (solid black line), is 
tracked closely by our numerical result for EDE (dashed red 
line). By contrast, the method of BDW (dot-dashed green line) 
predicts a si gnificant decrease in EDE cosmology. Our result 
for the iDoran fc Robbers! 1120061 ) EDE model (dotted blue line) 
has a s omewhat greater departure from the ACDM value than 
for the |Wettench| 1 2004) model, even at high redshifts. For the 
IDoran &: Robbers! 1 20061 ) model, the BDW approach gives 5 C — > 



(see Francis et al. 2008), so we do not plot that case. 

cosmologie:^]. By avoiding reference to turnover, we avoid 
this approximation. 

• By solving the exact equations, we can quantify the 
magnitude of any errors introduced by approximations that 
may lead to a simpler functional form for A(x). Without 
this solution we cannot properly test the validity of approx- 
imations. 

Results from our numerical calculation, compared to 
the method of BDW, are shown in Fig. [T] The EDE model 
is model (I) from BDW with parameters Q m = 0.325, 
wo = —0.93, and Q e — 2 x 10~ 4 . This is compared to a 
flat ACDM model with fi m = 0.3 as in B DW. We also com- 
pare t o different EDE model, p roposed bvlDoran fc Robbers! 
i|200d ). which, compared to the lWetterichl l|2004 ) model, has 
a greater EDE density at early times relative to its z ~ 2 
value. We show our results for this model, with £l m = 0.3, 
Wo = — 1 and fie = 0.05 (note that current data permits 
higher va lues of fi e in the IDoran fc Robbers! (120061) model 
than the IWetterichl (|2004l ) model, see IDoran et al.l (|2007l ) 
for details). The key result from Fig. [1] is that we find S c in 
EDE cosmology to be little changed compared to ACDM. 
This is significant, and in contrast to previous results. The 



1 From Eq. |5jl we see that the symmetry in the perturbation rise 
and fall times holds when the coefficient of y in the last term is 
time independent (then it gives the "development angle" solution 
y ~ cos 2 / 3 [C(t-t a )], where C is constant, cf. Peebles 1980). The 
coefficient is indeed time independent in three cases: 1) when A = 
(pure matter universe), 2) when g(x) is constant (which implies 
w(x) = —1: ACDM), and 3) when w(x) = —1/3 (equivalent to 
spatial curvature, i.e. a matter plus curvature universe such as 
OCDM). These were the cases most researchers were interested 
in until the late 1990s, and so symmetry has passed into the lore of 
linear collapse. However, it is not true for dark energy cosmologies 
with w 5^ — 1! 
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implications for the ST mass function are discussed in Sec- 
tion |3l The reasons for breakdown of the BDW solution are 
detailed in the Appendix. 

The dependence of S c on EDE is fit by 



trEDE I \ 

o c (ac) 



: A + [6(1 + w ) + c£l m - d] etc - efi e (10) 

where for the iDoran fc Robbers! (|2006l ) model, A = 1.6905, 
b = - 0.0183, c = 0.026 4, d = 0.0208 and e = 0.202, and 
for the IWetterichl l|2004h model, A = 1.6899, b = -0.0170, 
c = 0.0455, d = 0.0307 and e = 0.753. Both of these fits 
are good to ~ 0.1% in the range 0.2 < Q. m < 0.4, -1.2 < 
wq < -0.8, 0.1 < a c < 1.0 and < tt e < 0.05 for the 
IDoran fc Robbers! (I2006T ) model and < fi e < 1 x 10~ 3 for 
the lWetterichl (|2004|) model. By taking fl e = and w = — 1, 
Eg. PTOl returns the ACDM result. 



3 HALO MASS FUNCTIONS 

The results on the l inear collapse parameter solve a puzzle 
from iFrancis et alj (120081 ). where a marked difference be- 
tween the ST and J 01 mass functions i n ED E cosmologies 
was highlighted, and lGrossi fc Springell (|2008l ) where it was 
noted that if spherical collapse is ignored and a common 
value of 8 C = 1.689 is assumed instead the agreement be- 
tween these methods is restored. For the ST mass function, 
using spherical collapse arguments for S c , both studies relied 
upon the calculation from BDW. Using instead the method 
outlined in this study, we ha ve re-examine d the S T and J01 
mass functions. As found in lFrancis et alj (|2008l ). the EDE 
mass functions are not greatly altered compared to ACDM 
at z = 0, however the difference increases with redshift. 
The EDE mass functions as a ratio to ACDM at z = 1 for 
the same models as Fig. [T] are shown in Fig. [2] This result 
demonstrates that the basic agreement between the ST and 
the J01 mass functions is preserved, even when the spher- 
ical collapse motivated, cosmology dependent, S c is used. 
Thus not only the form, but also the conceptual basis of 
the ST mass function is valid for EDE. We note that at the 
high mass end (M > 10 13 Mq/K), the choice of calculating 
8 C or holding it fixed makes as big a difference as choosing 
betw een the ST and J01 for mulas. Neither the simula tions 
from lFrancis et all (|2008T ) nor lGrossi fc Springell (|2008) have 
sufficient accuracy at this high mass end to make any clear 
judgement about which choice better fits simulation data. 

The success of the J01 style mass functions, that are 
blind to the growth history of the universe (as opposed to 
the instantaneous growth factor), indicates that indeed the 
abundance of halos is insensitive to this. When the growth 
history is considered, via the alteration of 5 C in the ST mass 
functions, small, but not insignificant, differences between 
the ST and J01 predictions for the relative mass function 
in EDE and ACDM emerge. Future work with simulations 
containing sufficient volume to accurately probe the high 
mass range could discriminate between the two approaches 
to mass function fitting and determine whether the growth 
history affects halo abundances. 

While the abundance o f halos is unaffected by the 
growth history, IFrancis et al. I (|200Sl ) found that non-linear 
power at small scales (k > 1) is increased in EDE cosmolo- 
gies relative to ACDM. Since this part of the power spec- 
trum is dominated by the one-halo term in the halo model 
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Jenkins et al. (2001) 
Sheth-Tormen (1999) calculated 6 c 
Sheth-Tormen (1999) (5 C =1.6S9 




M (M Q /h) 

Figure 2. Comparison between the ST and J 1 mas s functions. 
The upper set o f curves are for the IWetterichl <2004) model and 
the lower set the lDoran fc Robb ers (2006) model. The solid black 



line shows the calculated ratio between the EDE and ACDM us- 
ing the J01 formula. The red dashed line shows the same ra- 
tio using the ST mass function for both models using our calcu- 
lated values of 8 C . The green dot-dashed line shows the ST ratio 
when <5 C = 1.689 i s use d instead of the calculated values. The 
IDoran fc Robber mass functions use linear power spec- 

tra, P(k), for the EDE and ACDM mode ls normalised to a com- 
mon erg today, while the lWetterichl (l2004l> mass functions use the 
models from BDW which are picked from a Monte-Carlo chain 
fitting current data and have some differences in normalisation 
and primordial spectral index. 



jCoorav fc Shethll2002l ). the internal density profile of halos 
in EDE will be different than in ACDM . This is also seen in 
the results of iGrossi fc Springell l|2008r ). 



4 DISCUSSION AND CONCLUSION 

In this letter we have demonstrated that the ST mass func- 
tion, when the corr ect 5 C is used, agrees with the J01 and 
IWarren et al mass functions. For reasonable param- 

eter values, 8 C in EDE models is not significantly altered 
compared to ACDM. This is in contrast to previous results 
for EDE, but agrees with a nalyses of other dyna mical dark 
energy models, for instance iMainini et al.l (|2003) . From Eq. 
(|3|, we can see that 5 C is defined by comparing linear to 
non-linear growth. If, in some cosmology, we find 8 C to be 
significantly altered compared to some other model, then 
this indicates that the difference between models must be 
altering the linear and non-linear growth differently. 

What kind of cosmology would have a significantly dif- 
ferent 8 C compared to ACDM? Some new physics must 
alter the linear and non-linear growth rates in different 
ways compared to ACDM. One key assumption we have 
made in our analysis is that dark energy is smoothly dis- 
tributed on the relevant length scales. This means that 
the dark energy density in the background universe and 
within the perturbation are not evolved independently. If 
dark energy perturbations were in fact important on small 
scales, or for instance if dark energy was non-minimally cou- 
pled to dark matter then the dark energy density within 
the perturbation will evolve differently to the background 
universe. In this case we would expect a more signifi- 
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cant alteration of 8 C , although careful determination of 
the linear growth as well as the non-linear spherical col- 
lapse taking any coupling or dark energy perturbations 
properly into account would be needed, see for instance 
llManera fc Motall2006l; lAbramo et al.ll2008l; iDutta fc Maori 



l2007l : lBasilakos fe VoriiJ2007l : lMainini fc Bonomettoll2006h 
As found in BDW, if S c is postulated to be signifi- 
cantly different to the ACDM value, then large differences in 
the abundance of collapsed objects would be seen, over and 
above any difference we might anticipate based upon mea- 
surements of the linear growth rate from the CMB, weak 
lensing and the large scale galaxy power spectrum normal- 
isation combined with knowledge of the expansion history 
from supernovae la data. Since some discrepancie s in struc- 
ture measurements may exist l|Fedeli et al.f l2008V it there- 
fore remains a question for further studies as to whether 
such observations potentially point to the non-linear growth 
in our universe not following the universal form and hint at 
added physics for dark energy or dark matter. 
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APPENDIX: THE BDW SOLUTION FOR A(x) 

The calculation of S c in BDW makes two key assumptions 
in order to approximate the early time evolution of the over- 
density in EDE cosmologies. In Eq. @, the second term is 
neglected on the assumption that it is small compared to the 
first, and the second order differential Eq. ([5]) is converted 
to an approximate first order equation. These approximate 
equations, for a flat universe, are 



[<j/x] 



1/2 



L0( (Q e LU 

V {1-U e )y 



1 1/2 



(11) 



(12) 



To examine the effects that these approximations introduce, 
we first determine ( using our numerical solution for A(oi), 
and then numerically integrate Eqs. (Ill|l and (I12|l in the 
early universe. If the approximations are good the results 
should match our exact calculation in the early universe. 

There is an important point that must be made using 
our approach. The boundary conditions of yi/yi = Xi/xi 
at at (the perturbation starts off co-moving with the Hub- 
ble flow) ensure that we have both a growing and decaying 
mode initially, such that S+ — (3/5)<5j where 5 = A — 1. 
At the initial time Xi, we have both growing and decaying 
modes and hence we need to apply the factor of 3/5 in or- 
der to extricate the growing mode only. The essence of the 
BDW approach is to approximate the early time trends of 
Eqs. @ and ©, which should return simply the growing 
mode. As expected from the results of Section [2] the inte- 
gration using Eqs. (|ll|l and (|12[) does not match our result, 
however it does reproduce the solution of BDW after the 
decaying mode has dissipated, indicating that these approx- 
imations cause the difference compared to our result (both 
approximations contribute comparable errors). The subse- 
quent manipulations of these equations in BDW are clever 
but the damage is already done. 

This paper has been typeset from a TgX/ DTjjX file prepared 
by the author. 



